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ABSTRACT 

We present an accurate analysis of the peculiar Horizontal Branch (HB) of the mas- 
sive Galactic globular cluster NGC 2808, based on high-resolution far-UV and optical 
images of the central region of the cluster obtained with HST. We confirm the mul- 
timodal distribution of stars along the HB: 4 sub-populations separated by gaps are 
distinguishable. The detailed comparison with suitable theoretical models showed that 
(i) it is not possible to reproduce the luminosity of the entire HB with a single helium 
abundance, while an appropriate modeling is possible for three HB groups by assuming 
different helium abundances in the range 0.24 < AY < 0.4 that are consistent with 
the multiple populations observed in the Main Sequence; (ii) canonical HB models are 
not able to properly match the observational properties of the stars populating the 
hottest end of the observed HB distribution, the so called "blue hook region." These 
objects are probably "hot-flashers," stars that peel off the red giant branch before 
reaching the tip and ignite helium at high T c g. Both of these conclusions are based on 
the luminosity of the HB in the optical and UV bands and do not depend on specific 
assumptions about mass loss. 
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1 INTRODUCTION 

Galactic globular clusters (GCs) are very old (age ~ 10-13 
Gyr) stellar systems populated by 10 5 -10 6 stars. Initially 
GCs were thought to be prime examples of simple stel- 
lar populations, i.e., populations of coeval stars born with 
the same initial chemical composition. That situation be- 
gan to change in the late 70's with the first spectroscopic 
surveys that disclosed star-to-star chemical variations of 
light element abundances within individual clusters (e.g., 
Cohen 1978). These variations are now known to be present 
in a large fraction of GCs and appear in the form of well de- 
fined anticorrelations between the abundances of C-N, Na- 
O and Mg-Al pairs, that cannot be explained in terms of 
evolutionary /mixing effects (see, e.g., Gratton et al. 2000, 
Carretta et al. 2008, 2009). In addition there is growing ob- 
servational evidence for helium (He) abundance variations 
within a few clusters which we discuss below. The presence 
of specific chemical patterns supports a scenario invoking 
multiple star formation events in some, perhaps most, GCs. 
These took place at the very beginning of the cluster evolu- 



tion on relatively short time-scales, of the order of 100 Myr. 
The younger generations could have been born out of gas en- 
riched by the winds of intermediate mass Asymptotic Giant 
Branch stars (AGBs; e.g., Ventura et al. 2009) or massive 
fast rotating stars (Decressin et al. 2007) formed during the 
first star formation episode. Dynamical simulations, e.g., by 
D'Ercole et al. (2008) and Decressin et al. (2008), have in- 
vestigated how these winds could have been retained by the 
clusters an incorporated into future generations of stars. 

One of the challenges facing self-enrichment scenarios in 
a few clusters is that they must be able to account for appre- 
ciable differences in the helium mass fraction (Y) between 
different generations, with Y reaching values up to twice the 
primordial abundance while at the same time not apprecia- 
bly altering the abundances of most other elements like iron. 
A notable exception is cj Centauri (NGC 5139), where the 
Main Sequence (MS) splitting detected by Bedin et al. 2004 
has been interpreted as due to He variation (with the richest 
population having Y ~ 0.38; Piotto et al. 2005) and where 
at least five distinct populations with different iron content 
and possibly different ages have been revealed (Pancino et 
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al. 2002, Ferraro et al. 2004, Sollima et al. 2005, Villanova 
et al. 2007, Calamida et al. 2009). However cjCen may not 
be a genuine GC, but rather the remnant of a dwarf galaxy 
partially disrupted by the interaction with the Milky Way 
(Bekki & Freeman 2003, Mackey & Van den Bergh 2005). 
A similar system (Terzan 5) harboring two distinct popula- 
tions with different iron content and (possibly) ages has been 
recently found in the Galactic Bulge (Ferraro et al. 2009). 
Similarly, this system is suspected to be not a genuine GC, 
but instead the remnant of a larger system that contributed 
to the formation of the Galactic Bulge. A milder evidence 
of internal spread of the metal content (A [Fe/H] ~ 0.1) has 
been observed also in M22 (Marino et al. 2009; Da Costa et 
al. 2009). 

In the context of the multi-population scenarios, one of 
the most puzzling cases is that of NGC 2808. The first pho- 
tometric observations (Harris 1974, Ferraro et al. 1990) re- 
vealed a Horizontal Branch (HB) with a very complex struc- 
ture. The cluster HB is well populated both at colours redder 
than the RR Lyrae instability strip and along its hot blue tail 
(BT) that covers a range of ~ 5 mag below the mean level of 
the instability strip. This morphology is not easily explained 
in terms of the cluster metallicity ([Fe/H] ~ —1-3), since it 
is well outside the common paradigm that links red HBs 
to metal-rich GCs and blue-HBs to metal-poor systems. In 
this sense NGC 2808 is similar to other very massive (and 
much more metal rich) GCs, like NGC 6388 and NGC 6441 
(Rich et al. 1997, Busso et al. 2007, Dalessandro et al. 2008). 
The cluster colour-magnitude-diagram (CMD) also shows 
puzzling discontinuities in the stellar distribution along the 
BT (Sosin et al. 1997, also Bedin et al. 2000, Castellani et 
al. 2006, Iannicola et al. 2009) similar to those found in sev- 
eral other GCs (Ferraro et al. 1998). 

D'Antona et al. (2005) first noticed a broadening of the 
NGC 2808 MS which was incompatible with photometric 
errors. They suggested that the MS consisted of two com- 
ponents with the same age and [Fe/H], but different initial 
He abundances, with Y ~ 0.4 for the He-rich component. 
The complex MS structure was confirmed by an accurate 
photometric and proper motion analysis with deep Hubble 
Space Telescope (HST) data (Piotto et al. 2007). They found 
that the MS of NGC 2808 splits into three sub-populations, 
all with age ~ 12.5 Gyr and Y ~ 0.248 for the red-MS, 
Y ~ 0.30 for the mean-MS and 0.35 < Y < 0.40 for the 
blue-MS. D'Antona et al. (2005) further hypothesized that 
the complex HB is connected to the MS components with 
different He abundances. Q However their HB morphology 
analysis was performed in the optical bands where the in- 
crease of the bolometric corrections at high effective temper- 
ature turns the HB into an almost vertical structure (rather 
than horizontal) when T > 10, 000 K. Along this vertical 
part of the HB, models with different initial He abundance 
tend to overlap, and the identification of sub-populations 
with different initial Y is not independent of assumptions 
about the amount of mass lost along the RGB phase. Ferraro 
et al. (1998) and Rood et al. (2008) suggest that the opti- 



1 Even though NGC 2808 is one of the most massive clusters and 
shows the presence of three distinct stellar populations charac- 
terized by quite different He abundances, ther is no evidence of 
appreciable [Fe/H] variations (Carretta et al. 2006). 



mal diagram for the study of blue HB stars is the (rriFisoBW , 
mFWOBW — rriF555w) CMD in the HST filter system. In this 
CMD the hottest HB stars are the most luminous and lie 
along almost horizontal sequences, whose luminosity is very 
sensitive to the initial Y abundance irrespective of the pre- 
cise value of the stellar mass. 

Here we present an accurate photometric study of the 
central regions of NGC 2808, based on high-resolution HST 
observations in far-UV and optical filters. With evolutionary 
models and constraints on the initial Y distribution coming 
from the analysis of the multiple MS, we investigate the 
complex structure of the HB of this massive GC. The pa- 
per is structured as follows. Section 2 describes the observa- 
tional data, their reduction and calibration, Sect. 3 presents 
an analysis of the observed properties of the cluster HB, 
and Sect. 4 details their interpretation in term of theoretical 
models. A summary and the conclusions follow in Sect. 5. 



2 OBSERVATIONS AND DATA REDUCTION 

We used a set of images (Prop. 6864, P.I. Fusi Pecci) cover- 
ing optical to far-UV wavelengths, obtained with the Wide 
Field Planetary Camera 2 (WFPC2) on board the HST. 
The WFPC2 is a mosaic made of four 800 x 800 pixels 
cameras, with angular resolutions of 0.046"/pixel for the 
Planetary Camera (PC) and ~ 0.l"/pixel for the three 
Wide Field Cameras (WF2, WF3 and WF4). The optical 
dataset consists of 5 x F555W images, three with expo- 
sure time t cxp — 100 sec and two with t cxp = 7 sec, plus 
4 x F336W images with exposure time t e x P = 1600 sec. 
The far-UV database comprises 3 x F160BW images with 
tcx P = 1200 sec each. All the images have the same point- 
ing, with the cluster centre roughly centred in the PC (see 
Figure 1). We combined the images using the IRAF task 
imcombine in order to improve the signal to noise ratio and 
decontaminate them from cosmic rays, which are particu- 
larly prominent in long exposure data. For a given passband 
and exposure time we adopted the resulting median frame 
as reference image for the data reduction. As done in our 
previous work (see for example Dalessandro et al. 2009), the 
data reduction of the optical images was performed using the 
ROMAFOT package (Buonanno et al. 1983), specifically de- 
veloped to perform accurate photometry in crowded regions 
(Buonanno & Iannicola 1984). Because of under-sampling 
problems related to the Point Spread Function (PSF) of 
stars in the far-UV images, we used aperture photometry 
performed with the publicly available software Source Ex- 
tractor (SExtractor, Bertin & Arnouts 1996). The adopted 
aperture radii were 2 pixels (corresponding to 0.2") for the 
WF chips and 3 pixels (~ 0.15") for the PC. For the WF 
cameras, because of strong vignetting problems affecting the 
upper right corners of all the F160BW images, the regions 
[550 < x < 800; 550 < y < 800] have been excluded from the 
analysis. In each filter for every star detected, the photomet- 
ric error has been defined as the standard deviation of all 
the measures obtained. 

The final star list includes all the sources detected in 
at least 2 filters. Hot stars detected in F160BW and not 
automatically found in the optical data were force-fitted in 
the optical images. This procedure added 10-15 HB stars 
in the F336W and F555W images. These stars are typically 
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Figure 1. Map of the WFPC2 data. The black cross corresponds 
to the cluster centre of gravity as determined in Sect. 12.11 The 
solid circle corresponds to the core radius (r c = 15.6"; Harris 
1996), while the dashed circle (r = 75") defines the area within 
which we limited our analysis. 

a magnitude fainter than the TO level in the optical CMD 
(Fig. 2). This approach allowed us to sample in a proper 
way both the cool sequences (like MS and RGB) thanks to 
the optical bands, and the very hot ones (extremely blue HB 
stars) because of the sensitivity of the far-UV bands to high 
temperatures. 

Our photometric strategy assured a very high degree of 
completeness in our samples. In particular, selecting stars 
where they are the brightest and then finding the corre- 
sponding stars in the plane where they are faint, guarantees 
an high completeness level in each plane. The F160BW im- 
ages are not affected at all by crowding problems, allowing a 
quite trivial detection of all hotter objects. As a consistency 
check we counted the intermediate-temperature HB samples 
in both the optical and UV planes. The samples are virtu- 
ally identical (240 stars in the optical and 243 in the UV 
catalogue) . 

2.1 Photometric calibration, astrometry and 
centre of gravity 

The WFPC2 photometric catalogue has been placed onto 
the absolute GSCII astrometric system by using the pro- 
cedures described in Ferraro et al. (2001, 2003). To max- 
imize the number of stars in the WFPC2 Field of View 
(FOV) adopted as reference for the roto-traslation proce- 
dures, we used the photometric catalogue described by Fab- 
bri et al. (2010) that was obtained with the Wide Field 
Imager (WFI) at the 2.2 m ESO-MPI La Silla telescope, 
as secondary astrometric standard list. With the help of 
the CataXcorr package, developed by Paolo Montegriffo at 
the Bologna Astronomical Observatory, we found some hun- 
dreds of stars in common between the WFPC2 and the WFI 



Figure 2. Optical (m_F555w , m F336w — m F555w ) CMD. We also 
show the boxes used to select RHB and BHB stars. 

catalogues. These were used to obtain a very accurate as- 
trometric solution. At the end of the procedure the typical 
error is ~ 0.2" both in Right Ascension (a) and Declination 
(6). 

All the magnitudes were transformed to the VEGA- 
MAG photometric system using the procedure described in 
detail by Holtzman et al. (1995) with the gain settings and 
ZeroPoints listed in Tab. 28.1 of the HST data handbook. 
The resulting CMDs are shown in Figs. 2 and 3. 

We determined the centre of gravity (C gra v) of 
NGC 2808 with an iterative procedure (Montegriffo et 
al. 1995) by averaging the absolute positions (a and S) of 
stars lying within r = 10" from a given point, starting from 
the centre reported by Harris (1996). As done in previous 
papers (see Lanzoni et al 2007 for an example) we selected 
four different sub-samples of stars with different magnitude 
limits (rriF555w = 20, 19.5, 19 and 18.5) to avoid incomplete- 
ness effects in the central region. The average of the four 
derived values has been adopted as C gra v and it is located 
at ajaooo = 9 h : 12 m : 03.06 s and <5 J20 oo = -64°51'48.82", 
~ 3.5" northeast from the centre quoted by Harris (1996). 



3 THE OBSERVED HORIZONTAL BRANCH 
POPULATIONS 

In the canonical view, the HB is expected to be populated 
by stars with different masses, decreasing from cold to hot 
effective temperatures, all undergoing central He-burning 
and shell H-burning. Given that the mass of the He-core 
at the start of the HB phase (Zero Age Horizontal Branch- 
ZAHB) is almost constant (M c ~ 0.5M©) and equal to the 
core mass at the He-flash, the location of each star along 
the HB is essentially determined by the mass of the sur- 
rounding envelope that survives the mass-loss episodes dur- 
ing the RGB phase. From the observational point of view, 
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Figure 3. Bottom panel. Far-UV (»hfi60-BW > m Fi60BW ~ 
m F555w) f° r M80 which we adopt as a "normal" blue tail clus- 
ter. Data are from Ferraro et al. (1998). The solid grey line is a 
fiducial line that reproduces the lower envelope of the HB star 
distribution. The colour and magnitudes have been shifted to 
compensate for the difference in reddening and distance mod- 
ulus between M80 and NGC 2808 as determined by demand- 
ing that lower envelopes of the two BHBs match in the interval 
— 1 < (rap iaoBW ~ m F555w) < 2. Middle-panel. The far-UV 
CMD for NGC 2808. Different symbols denote different HB sub- 
populations: open circles are BHB stars, open triangles are EHB 
stars and open pentagons are candidate BHk stars. The dashed 
vertical line marks the hot edge of the distribution of HB stars 
in M80. Upper upper panel. Far-UV star count distribution of 
NGC 2808. 



the first physical parameter affecting the HB morphology 
is the metallicity. Metal rich clusters usually have red HBs, 
while the metal-poor ones tend to have HBs populated at 
bluer colours. However, at least one additional, or second, 
parameter is required to explain the observed colour distri- 
butions of HB stars in Galactic GCs. This second parameter 
problem has received much attention during the last four 
decades (see for example Gratton et al. 2010 and references 
therein). Age has been the most commonly discussed second 
parameter (Lee et al. 1994; Dotter et al. 2010; Gratton et al. 
2010). Very recent work by Gratton et al. (2010) and Dotter 
et al. (2010) shows that age as a global second parameter can 
account for general HB parameters such as mean colour and 
length of the HB. However, HB morphologies show a diver- 
sity which cannot be fully characterized by mean properties. 
Sets of clusters with similar metallicities (M3/M13/M80, see 
Ferraro et al. 1998; 47 Tuc/NGC 6388/NGC 6441, Rich et 
al. 1997) have HBs distinguished by gaps and long blue tails. 
Features like these have prompted discussion of additional 
parameters like mass-loss efficiency during the RGB (Rood 
et al. 1993; Ferraro et al. 1998; Dotter 2008) or initial He 
abundance differences (D'Antona et al. 2002). The situation 
is still more complicated because deviations from the gen- 
eral behaviour arise due to internal second or third param- 



eters that may vary from cluster to cluster. In this context 
NGC2808 is a particulary interesting case which certainly 
requires an accurate and detailed comparison between ob- 
servations and state-of-the-art theoretical models. As shown 
in Figs [2] and the HB of NGC 2808 is characterized by 
a well populated red clump and, on the blue side, a well 
populated and extended blue sequence reaching down to 
1TIF555W ~ 21.5, almost 2 magnitudes fainter than the MS 
turn-off (TO). 

By combining observations in optical and far-UV bands 
we can sample HB stars according to their temperatures and 
spectral type. The strategy adopted here is to select stars 
along the HB in the most appropriate CMD according to 
their temperature. Hence, the red clump stars have been se- 
lected in the optical CMD (mFsssw, rn,F336w — thfsssw), 
that is more sensitive to relatively cool stars, while the 
hottest HB stars were selected in the far-UV CMD, where 
they are the most luminous sources. These selection criteria 
allow also a more accurate distinction between pure HB and 
post-HB stars. The intermediate-temperature HB popula- 
tion has been selected in both CMDs, and it was used a link 
between the cool and the hot part of the HB, as well as to 
define normalization criteria and the relative incompleteness 
of the two samples. Two under-populated narrow regions (or 
"gaps") are apparent in Figures 2 and 3. These gaps were 
previously identified by Bedin et al. (2000; hereafter BOO) 
at V = 18.6 and V = 20 and the location in our CMD (at 
rriF555w = 18.55 and rriF555w ~ 20.05) is well consistent 
with previous estimates. The two gaps split the HB in 4 
subgroups (if we also consider the standard subdivision at 
the cool and hot sides of the RR Lyrae instability strip) , that 
were named RHB, EBT1, EBT2 and EBT3 by BOO (see also 
Castellani et al. 2006; Iannicola et al. 2009, hereafter 109). 
For sake of homogeneity with previous work of our group, 
we will use the nomenclature introduced by Dalessandro et 
al. (2008), and we will denote the four sub-populations as 
red HB (RHB), blue HB (HB), extreme HB (EHB) and blue- 
hook (BHk), respectively. The physical properties of these 4 
sub-groups will be discussed in the following section. 

Because of well known vignetting problems affecting the 
upper-right corner of all F160BW images (see Sect. [2}, we 
restrict our analysis to stars with distances r < 75" from 
Cgrav The selected populations are shown with different 
symbols in Figs. 2 and 3. The total HB sample in our cat- 
alogue is 616 stars: 256 RHB, 243 BHB, 66 EHB and 51 
BHk sources. They correspond, respectively, to (41 ± 3)%, 
(39 ± 3)%, (11 ± 2)% and (9 ± 1)% of the total. These ratios 
are in agreement with those reported by 109 and are com- 
patible with those quoted by BOO. We found 4 RR Lyrae in 
common between the variable star list published by Samus 
et al. (2009) and Corwin et al. (2004) and our BHB sam- 
ple. Given the small number, their inclusion doesn't effect 
the comparison with theoretical models performed in the 
following sections, nor the results obtained in this work. 

To fully appreciate the oddity of NGC 2808's HB, it 
is useful to compare it to a "normal" blue tail cluster, 
M80. Of the clusters we've studied in the far-UV, M80 
has the most extended BT. Its CMD along with that of 
NGC 2808 is shown in Figure [3] To facilitate comparison 
the colour and magnitude have been shifted to compen- 
sate for the difference in reddening and distance modulus 
between M80 and NGC 2808 as determined by demanding 
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Figure 4. Cumulative radial distribution of the four different HB 
sub-populations. 

that lower envelopes of the two BHBs match in the interval 
— 1 < (m,Fi60BW — rriF555w) < 2 and added a fiducial. The 
fiducial marks the lower envelope of the M80 HB0 There 
are two dramatic differences: first, the HB of NGC 2808 is 
more extended toward high temperatures than that of M80; 
second, if the fiducial is linearly extrapolated to higher tem- 
perature many of the NGC 2808 stars lie below the fiducial. 
To fit the observed NGC 2808 sequence the fiducial would 
have to hook downward. While these differences become evi- 
dent in far-UV CMDs, they are barely recognizable in purely 
optical diagrams. 

The spatial distribution of a stellar population can give 
information about its origin. Figure [4] shows the cumula- 
tive radial distributions of the 4 HB subgroups. The RHB, 
BHB, EHB and BHk stars have similar radial distributions 
within the WFPC2's FOV. A Kolmogorov-Smirnov (KS) 
test shows that there are no statistically significant dif- 
ferences between the radial distributions of these 4 sub- 
populations. The strongest discrepancy is found between 
RHB and BHB stars, but it is significant only at the ~ 1.3<r 
level. 109 and BOO found similar results over the area covered 
by our data. 109 did find significant changes in the ratios of 
the different HB subgroups between the centre and their 
outer zone (r > 90"). 



4 THEORETICAL ANALYSIS 

We have compared our photometry with an appropriate set 
of BaSTI a-enhanced HB tracks and isochrones (Pietrin- 
ferni et al. 2006) for [Fe/H] = -1.31 close to the old Zinn 

2 Using ZAHBs of different metallicity we have determined that 
it is reasonable to use the shifted M80 fiducial for NGC 2808 
despite the difference in metallicity. 
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Figure 5. (m P555w , m P336w - m F555w ) and {m F336w , 
m F336W ~ m F555w) CMDs for the cluster RHB population, and a 
portion of the cluster RGB. The box in the lower panel marks the 
region used in the comparison with synthetic HB diagrams (see 
text for details). The solid lines are the theoretical ZAHB and 
the isochrones at t = 12Gyr and [Fe/H]= -1.31 with Y = 0.248, 
shifted by the distance and extinction determined for the cluster. 



& West (1984) value of —1.36. Current spectroscopic es- 
timates are larger ([Fe/H] = —1.18, Rutledge et al. 1997; 
[Fe/H] = -1.10 or -1.15, Carretta et al. 2009). A change of 
0.2 in [Fe/H] shifts a ZAHB by 0.08 mag. A comparison with 
the observations can be made with a compensating shift in 
distance modulus. 

We considered the entire set of initial He abundances 
available in the BaSTI database, i.e. Y=0.248, 0.300, 0.350 
and 0.400. 

The bolometric corrections to the WFPC2 bands have 
been determined using the a-enhanced spectral library by 
Castelli & Kurucz (2004) — the same one adopted by Pietrin- 
ferni et al. (2006) — and following the procedure outlined by 
Girardi et al. (2002). This method allows a straightforward 
calculation of the effect of extinction in the adopted fil- 
ters, and its dependence on both effective temperature (T e ff) 
and surface gravity. We adopted the Cardelli, Clayton & 
Mathis (1989) extinction law, with Ry = 3.1. 

The spectroscopic study by Pace et al. (2006) showed 
that in NGC 2808 HB stars hotter than ~ 12000 K have sur- 
face [Fe/H], [Cr/H] and [Ti/H] values around solar or even 
higher (up to [X/H] ~ 1, where X is the abundance of one of 
these three metals) due to radiative levitation in the atmo- 
sphere. In the ZAHB models and HB tracks here adopted we 
mimic this effect by applying bolometric corrections appro- 
priate for [Fe/H] = 0.0 when T e g is between 12,000 K and 
13,000 K, and, when the temperature is higher, bolometric 
corrections appropriate for [Fe/H] = 0.5, based on the spec- 
troscopical evidence provided by Pace et al. (2006, see their 
Figure 4). 

This is in principle a crude approximation, but extended 
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Figure 6. Comparison between observed (grey histograms) and 
theoretical (dashed histograms) star counts for the RHB popula- 
tion, as a function of mpszew and rriFBBBW magnitudes. 

(in both mass and initial chemical composition) grids of HB 
stellar evolution and atmosphere models that include consis- 
tently the effect of radiative levitation (see , e.g. Hui-Bon- 
Boa, LeBlanc & Hauschildt 2000, and Michaud, Richer & 
Richard. 2008 for first results on this topic) are still lacking. 
Also, it appears that additional mixing processes (e.g. merid- 
ional circulation) must be coupled to the effect of radiative 
levitation, in order to explain the abrupt disappearance of 
the surface abundance anomalies at T e //11000-12000 K (see, 
e.g., Quievy et al. 2009). 

It is clear that our way to simulate the effect of radiative 
levitation is just an approximation, that cannot take into 
account a possible stratification of the metal abundances in 
the atmosphere and below the photosphere. In the follow- 
ing we will however show how this treatment seems to be 
reasonably appropriate. 

4.1 Determination of the cluster reddening and 
distance modulus 

As a first step we determined a reference distance modu- 
lus and extinction by fitting our HB models to the RHB. 
First, we produced synthetic samples of HB stars by follow- 
ing the techniques pioneered by Rood (1973; see also Ferraro 
et al. 1999) and compared separately the magnitude distri- 
butions of synthetic RHB stars in the F336W and F555W 
filters to the observed ones for the stars populating the box 
displayed in the lower panel of Figure [5] RHB stars are too 
faint in the F160BW filter to provide any additional con- 
straint. 

We used the set of models with Y = 0.248 for the RHB 
stars. The synthetic samples have four free parameters: ex- 
tinction, distance modulus (m — M)o, mean HB stellar mass 
(Mhb) and its dispersion C(jvf HB >- We started by randomly 
selecting the stellar mass Mhb from a Gaussian distribution 




-2.0 -1.0 0.0 1.0 



1U F336W U1 F555W 

Figure 7. Optical CMD of the cluster HB compared to theoret- 
ical ZAHBs for Y=0.248, Y=0.30 and Y=0.40 (see labels). The 
dashed lines show the effect of neglecting radiative levitation in 
the ZAHB bolometric corrections when T fj is larger than 12000 K 
(see text for details). The typical behaviour of the photometric 
errors is also shown as a function of magnitude. 



centred around a given guess value of (Mhb), with la dis- 
persion o"(m H b>- The WFPC2 magnitudes of the synthetic 
star were determined according to its position along the 
appropriate HB track with mass Mhb (interpolated, when 
necessary, among the available set of tracks in the BaSTI 
database) after an evolutionary time t. We determined t as- 
suming that the stars reach the ZAHB with a constant rate. 
We employed a flat probability distribution ranging from 
zero to £hb , where £hb is the time spent along the HB (we 
consider as the end of the HB phase the time when the cen- 
tral He-abundance drops to zero). The star with the lowest 
mass has the longest central He-burning lifetime (Castellani 
et al. 1994, Cassisi et al. 2003) and sets £hb ■ This implies 
that for some masses the randomly selected value of t will 
be longer that their HB lifetime, or in other words that they 
have already evolved to the following evolutionary stages. 

We modifed these synthetic magnitudes accounting for 
distance modulus (m — M)q and extinction. We used as 
broad guidelines the results from previous investigations 
(e.g. Bedin et al. 2000), exploring a range of extintion values 
E(B - V) = 0.15-0.20, and distance modulus (m - M) = 
15.0-15.5. We finally added a Gaussian random error with 
la dispersion, equal to that obtained from the data reduc- 
tion process. We then compared the resulting synthetic dis- 
tribution of m,F555w and m,F33ew values to the observed 
one, and modified the values of the four free parameters un- 
til the mean magnitude and la dispersion equal the observed 
ones in both mFBBBW and mFsmw- 

For each choice of the four free parameters we have 
produced 200 synthetic distributions, each of them with the 
same number of objects as in the observational sample. The 
synthetic star counts as a function of magnitude displayed in 
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Figure 8. As in Fig. [7] but for the (m F160BW , m Fle0BW — 
m F555w) CMD. The RHB population is not displayed because 
too faint in the F160BW filter. The dashed-dotted line represents 
the hot-flasher model used to reproduce the BHk stars, as de- 
scribed in Section 14.31 The behaviour of the photometric errors 
as a function of colour is also shown. 

the rest of the paper are average values from these ensembles 
of simulations. 

We found Afsbsw = 0.52 (corresponding to E(B — V) = 
0.17), (m - M) = 15.23, (Mhb) = 0.69M© and ct (Mhb) = 
O.O15M . 

The theoretical ZAHB and the RGB of a 12 Gyr 
isochrone with Y — 0.248 shifted according to the best-fit 
distance modulus and extinction are shown in Figure [5] Fig- 
ure[5]shows the magnitude distributions for the RHB stars in 
the observed (shaded histogram) and synthetic (dashed his- 
togram), for the best fit parameters obtained above. We also 
checked — comparing theoretical and observed star counts 
along the RGB — that the magnitude of the so-called RGB 
bump is well reproduced in both rriFsssw and rriF33ew by 
the 12 Gyr theoretical isochrone, for the combination of dis- 
tance modulus and extinction determined from the HB. 

By considering the initial mass of the objects populating 
the tip of the RGB in the 12 Gyr isochrone, we found that 
the mean amount of mass lost along the RGB by the Y = 
0.248 population is (AM) = 0.15 Mq. This value is very 
similar to the results of the synthetic modeling of the RHB 
stars by D'Antona & Caloi (2004), who obtained (AM) = 
0.13 Mq and a dispersion &(m hb ) = 0.015 Mq, using their 
scaled solar isochrones with approximately the same [Fe/H] 
of our Q-enhanced models (and slightly lower Y). 

4.2 Multiple-He populations along the HB 

To investigate the nature of the BHB, EHB and BHk popu- 
lations we used the most appropriate CMDs: the (mF33ew, 
mF336w — rnF555\v) CMD has been used as a bridge be- 
tween the RHB and BHB populations, and the (rriFimBW , 



TriFWOBW — mF555w) CMD to characterize the BHB, EHB 
and BHk populations. Prompted by the multiple MS de- 
tected in the cluster (Piotto et al. 2007) we used models with 
different He content for various components of the HB. Fig- 
ure[7]shows the (ttif336W, itif336W — wifsssw) CMD of RHB 
and BHB stars. Using the distance modulus and extinction 
obtained in the previous section, we also display three ZA- 
HBs with different Helium abundances Y = 0.248, 0.30 and 
0.40 respectively. For each ZAHB we show two sequences, 
one that mimics the inclusion of levitation as described 
above (solid line), and one without levitation (dashed line). 
We see that: 

(i) Applying metal-rich bolometric corrections to mimic 
radiative levitation significantly affects the location of the 
ZAHB in the high temperature regime, i.e. for (rn,F336w — 
w.F555w) < 0. In particular, the "levitation" ZAHB is sys- 
tematically bluer (S(mF336W — iriF555w) ~ 0.4) than the 
"no-levitation" ZAHB (or the ZAHB level is brighter at fixed 
colours) . 

(ii) The BHB stars define a sort of "knee" in the stel- 
lar distribution and are significantly brighter than the RHB 
clump in this plane. Figure [7] shows very clearly that even in- 
cluding radiative levitation, the magnitude level of the BHB 
sub-population is not reproduced by the ZAHB with canon- 
ical helium abundance (Y = 0.248): the BHB appears sig- 
nificantly (0.3-0.4 mag) brighter than the canonical ZAHB. 
On the other hand, the "levitation" ZAHB with Y = 0.30 
matches very nicely the lower envelope of the BHB star dis- 
tribution. 

(iii) The cooler BHB stars, (mF336w — wifsssw) > 0, 
are not affected by levitation, yet they still lie significantly 
above the canonical Y = 0.245 ZAHB. They are well fit 
by the Y = 0.30 ZAHB. The conclusion that Y for the 
BHB is ~ 0.05 greater than Y of the RHB is independent 
of levitation and our treatment thereof. 

The situation is less clear for EHB and BHk stars 
because their distribution becomes almost vertical in this 
CMD, and photometric errors increase significantly the scat- 
ter at faint magnitudes. For these reasons we used the 
(rriFieoBW, m F ieoBW - rn F 555w) CMD in order to study 
the properties of these sub-populations. Figure [5] shows the 
CMD of the three bluest sub-populations, together with the 
ZAHB for different initial values of Y, with and without 
inclusion of levitation. It is interesting to note that in the 
F160BW filter the inclusion of levitation goes in the opposite 
direction than in optical filters, making the ZAHB fainter at 
a given colour. Additional interesting features emerge from 
this figure: 

(i) Just as it was in the (m F 336W, m F 336w - m F 555w) 
CMD, the BHB population is matched by a ZAHB with 

Y — 0.30 including the effect of levitation. 

(ii) The EHB population is consistent with a ZAHB with 

Y — 0.40 including levitation. 

(iii) The BHk stars are too hot to be reproduced by any 
ZAHB models, hence they require a more separate discussion 
(see Section 4.3) 

From this comparison we can safely conclude that 
three (out of four) of the sub-populations along the HB of 
NGC 2808 can be identified as the progeny of three dis- 
tinct stellar populations with different initial helium abun- 
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dances in the range Y = 0.24-0.40. The specific values used 
(Y = 0.248, 0.30, 0.40) were set by the model data base, and 
we have not attempted to determine values of Y which best 
fit the data. The adopted values are consistent with what 
found from the photometry of the MS (Piotto et al. 2007). 
This conclusion is independent of any assumption about the 
mass-loss efficiency of their RGB progenitors. Our findings 
are also in agreement with He abundances reported on the 
basis of different indicators along the RGB (like difference 
in colour or T e a) by Bragaglia et al. (2010). The consistency 
with these results also justifies our choice (in principle ar- 
bitrary) to fix the distance (and extinction) by matching 
the Y = 0.248 HB population to the RHB. Had we tried to 
match these "normal-He" models to, e.g. the BHB popula- 
tion, it would have been impossible then to match the RHB 
stars with models with any reasonable value of Y. 

As a second step of our analysis, we determined the 
mean mass of the BHB and EHB sub-populations, by pro- 
ducing synthetic HBs in the same way as described be- 
fore, and matching the mode of both the m,FS36W and the 
(mFWOBW — mF555w) distributions. The synthetic sample 
compared to the BHB component has been calculated using 
the Y = 0.30 HB tracks, while for the EHB we employed 
the Y = 0.40 tracks, including the effect of levitation in 
both cases. Figure [9] compares the observed (shaded his- 
togram) and theoretical (solid histogram) star counts as a 
function of rriF336w and the (rriFieoBW — itifsssw) for the 
BHB and the bluer EHB component. An indirect estimate 
of mass-loss experienced by each sub-population can then be 
obtained by comparing the mean mass on the HB and the 
mass at the TO (for an age of 12Gyr). A summary of the 
values obtained for RHB, BHB and EHB sub-populations 
is reported in Table 1. Under these assumptions we found a 
mean amount of RGB mass-loss ranging from (AM) = 0.15 
for both the RHB and EHB stars, to (AM) = 0.20 for the 
EHB stars. 

Finally, in order to verify whether the RHB, BHB and 
EHB populations are the progeny of the red, mean and blue 
MS, respectively, we have compared the star counts along 
the HB, with those reported by Piotto et al. (2007) along the 
MS. In our sample we found: (A?RHB/iVBHB)obs = 1-1 ± 0.1, 
and (AfRHB/^VEHB)obs = 3.9±0.7. Starting from the observed 
MS number counts, and after applying corrections account- 
ing for the different evolutionary lifetimes along the RGB 
and among the HB populations (i.e., how the lifetime varies 
as a function of Y and mass) the expected number of stars 
along the HB would be: (./Vrhb/-/Vbhb)cxp = 3.7 ± 0.4 and 
(■ATrhb/ATehb^xp — 3.3 ± 0.5. Clearly, while the two values 
of (A^rhb/A^ehb) are compatible, it is not possible to recon- 
cile the observed and the expected values of (A^rhb/A'bhb). 
The ACS data used by Piotto et al. (2007) are strongly off- 
centred (r ~ 200") with respect to C grav . As shown by 109, 
the ratio -/Vrhb/A^bhb roughly doubles for r > 90". This 
may be a possible explanation for the observed discrepancy. 
However the progenitor masses for RHB and BHB are re- 
spectively 0.84 Mq and 0.76 M (see Table 1). Dynamical 
evolution of the cluster should cause the A^rhb /Abhb to de- 
crease as r increases, just the opposite of what was found by 
109. A more detailed study of radial variations in HB sub- 
populations extending over the full cluster could be illumi- 
nating. It would also be quite interesting to check whether 
the triple MS varies with position in the cluster. 
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Figure 9. Comparison between observed (grey histograms) and 
theoretical (solid histograms) star counts for the BHB (matched 
with models with y=0.30) and EHB populations (matched with 
models with y=0.40) as a function of the m.F336VK magnitudes 
(upper panel) and {mpimBW — m F555w) colour (lower panel). 



4.3 The Blue Hook population 

The CMDs in Figs. [7] and [8] show very clearly that BHk 
stars are not reproduced by any of the theoretical mod- 
els, this sub-population is systematically hotter than the 
hottest point along the ZAHBs. To shed more light on the 
nature of these stars, in Fig. [TO] we compare the observed 
number count distributions for the BHk sub-population in 
the F160BW, F336W and F555W magnitudes, with a 
synthetic sample with Y = 0.40, a mean mass (A/hb) = 
0.461Mq (i.e. the minimum possible mass of HB objects 
for our models at this metallicity) and zero dispersion. It is 
impossible to match the observed distributions in all three 
filters. 

This inability of HB models to reproduce the BHk pop- 
ulation suggests that these are not genuine HB stars, and it 
leads us to identify these objects as hot-flashers (D'Cruz et 
al. 1996; Brown et al. 2001, 2010; Cassisi et al. 2003; Moehler 
et al. 2004). 

Hot-flashers are stars that experience strong mass-loss 
during the RGB phase, leave the branch before the occur- 
rence of the He-flash, and move quickly to the He-core white 
dwarf cooling curve, where they experience a He-flash under 
conditions of strong electron degeneracy in their core. This 
scenario was envisaged earlier by Castellani & Castellani 
(1993) who named these peculiar stellar objects "red giant 
stragglers." Stars that ignite helium on the white dwarf cool- 
ing sequence have a much less efficient H-burning shell, and 
will undergo extensive mixing between the He-core and the 
H-rich envelope (Cassisi et al. 2003). This makes the progeny 
of hot-flashers hotter and fainter than genuine HB stars (i.e. 
stars that ignite He along the RGB), describing a kind of 
hook in UV CMDs that gives the name to this class of ob- 
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Figure 10. Comparison between observed (grey histograms) and 
theoretical (dashed histograms - models with y=0.40) star counts 
for the BHk population. The synthetic sample reaches the bluest 
possible colours allowed by the theoretical models. 

jects. The presence of such peculiar sources in NGC 2808 
has been already noted by pure UV analysis by Brown et 
al. (2001) and Dieball et al. (2005). As shown in Figure H 
the location of BHk stars is indeed reproduced by a hot- 
flasher model (dashed-dotted line) from Cassisi et al. (2003) 
for a mass M = O.489M 

The use of hot-flasher models with different total mass 
would not greatly modify the outlined scenario since all 
hot-flasher models tend to cluster in a very narrow re- 
gion of the CMD, regardless of their total mass (Brown 
et al. 2001; see also Miller Bertolami et al. 2008). Also, 
the effect of the initial chemical composition on the lo- 
cation in Te« is minor, at least at subsolar metallici- 
ties (see Miller Bertolami et al. 2008). As a consequence 
it is not trivial to understand which Y sub-population 
produces the BHks. Potentially star counts derived from 
MS analysis might shed light on this point. Hypoth- 
esizing three different progenitors for BHks we obtain 

(A" R HB+BHk/A r BHB)obs = 1.3 ±0.1, (-/VrHb/ ^VaHB+BHlcjobs = 

0.9 ±0.1 and (A r RH B/A r EHB+BHk)ob s = 2.2 ± 0.4. In all cases 
adding the BHk stars to a given HB sub-populations makes 
the disagreement with the ratios expected from Piotto et al. 
(2007) worse. 

3 The bolometric corrections used in the hot-flasher model dis- 
played in Figure \E\ have been obtained from scaled-solar model 
atmospheres, accounting for the effect of levitation as described 
in Section. 4. As shown by Brown et al. (2001), scaled solar spec- 
tra do not accurately represent the peculiar atmospheres of hot- 
flashers, which are expected to have strongly enhanced He and 
C abundances. Therefore, for the hot-flasher model one should in 
principle use more appropriate — but not yet publicly available — 
bolometric corrections, that however, should not alter the basic 
qualitative conclusions of our analysis. 
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Figure 11. The observed distribution of T eS for NGC 2808 HB 
stars. Grey histograms are the temperature distribution for ob- 
served stars obtained as discussed in Section 14.41 Dashed black 
lines rapresent the temperature distribution of simulated stars. 
The four well separated groups of stars are matched with models 
with, respectively - moving from the lowest to the highest T e yy 
- y=0.248, 0.30, 0.40, and the adopted track of hot-flashers. 

Bailyn (1995) has suggested that hot HB stars in a num- 
ber of clusters have a binary origin. On the other hand a 
recent spectroscopic analysis (Moni Bidin et al. 2006) of a 
sample of hot HB stars in NGC 6752 found no close binaries. 
Neither of these studies have direct application to the case 
of NGC 2808. Bailyn's arguments apply to sdB stars which 
are probably EHB stars (Dorman et al. 1993). Our prelim- 
inary results indicate that NGC 6752 has no BHk stars, so 
the Moni Bidin result also applies to EHB stars. Our results 
here clearly show that EHB and BHk stars are not the same 
thing. One must be quite precise in defining these hot pop- 
ulations. Specifically in NGC 2808, Piotto et al. 2007 tenta- 
tively suggest that their observed binary sequence might be 
connected with the BHk stars, and they note that the bi- 
nary fraction and BHk fraction are similar. Most GCs have 
binary stars, and most of these contain exotic populations 
like BSS which are produced by binary evolution. On the 
other hand, BHk stars are found in only a very few GCs. If 
binaries are producing BHk stars, they should be found in 
most GCs. So far genuine BHk stars have been found only 
in the most massive clusters most of which have multiple 
populations. This must be telling us something about their 
origin. 

4.4 The HB temperature distribution 

The combined use of far-UV and optical filters allows us 
to properly visualize the different HB sub-populations and 
reveal the existence of gaps, avoiding artificial features gen- 
erated by the saturation of the chosen colour index as it 
becomes insensitive to changes in T e g. The bottom panel 



10 E. Dalessandro et al. 



of Figure [9] shows, for example, how the HB distribution in 
the (rriFieoBW — "Ifsssvk) colour separates very clearly the 
EHB and BHB sub-populations. In order to obtain the dis- 
tribution in temperature of the different HB groups, we first 
derived T e ff from the observed and dereddened colours by 
interpolating a spline along the colour- T e // relation given 
by the appropriate theoretical ZAHBs. The result is shown 
as grey histogram in Figure [TT1 

This approach neglects the post-ZAHB evolutionary effects 
on the colour-T e / / relation which however are small. In fact 
the derived distributions have been found to be fully com- 
patible with those directly obtained from the synthetic mod- 
els (dashed histograms in Figure [TTj) . 

The T c ff distribution displayed in Figure [TT1 appears clearly 
multimodal, as expected. Each peak observed in this dis- 
tribution belongs unequivocally to only one of the 4 sub- 
populations. The effective temperatures spanned by the 
HB stars in NGC 2808 cover the range 5,000 < T cff < 
50, 000 K. The RHB has a distribution peaked at Trhb ~ 
5,400K, the BHB peaks at T BH b ~ 12, 600 K, the EHB 
at Tehb ~ 24, 000K, while BHk stars are centreed around 
T BH k ~ 40, 000 K. The gaps Gl and G2, that we have found 
at m,F555w = 18.55 and tofsssw = 20.05, are located at 
T G i ~ 20, 000 K and T G2 ~ 31600 K, respectively. The tem- 
peratures of these two gaps (Tg2 in particular) are slightly 
hotter than those (~ 17, 000 K and ~ 25, 000 K) estimated 
by Sosin et al. (1997). These latter values are based on mod- 
els with a constant and standard (Y = 0.236) He abun- 
dance. Moreover, since they were unable to fit both the op- 
tical and near UV data with the same distance modulus 
and reddening values, they applied colour shifts (for example 
A(B - V) = 0.07 towards the red in the optical CMD) that 
may explain the discrepancies in the gaps' temperatures. 
Our estimate for Tgi and Tg2 are much more compatible 
with the temperatures estimated by the spectroscopic ob- 
servations of a sample of EHB and BHk stars by Moehler 
et al. (2004; see their Table 4 and their Figure 2). Their 
BHk candidates have typically T e g > 35, 000 K, while the 
few EHB stars observed are distributed in the temperature 
interval 20, 000 < T c « < 30, 000K . Our estimated highest 
temperature along the HB is substantially higher than the 
value (~ 37, 000K) reported by Recio-Blanco et al. (2006) 
from optical CMDs. This once more demonstrates the diffi- 
culties of deriving the properties of extremely hot HB popu- 
lations from pure optical CMDs. In the optical CMD a star 
moves down the BT as the bolometric correction increases 
with Tgfs or the stellar luminosity decreases. For BHk stars 
both are changing and optical colours cannot give a reliable 
T c ff . However it is important to note also that the tempera- 
tures estimated for BHks are affected by large uncertainties 
that are not easily quantifiable. They are due to the lack 
of a proper treatment of the peculiar chemical composition 
of their atmospheres and to the photometric errors that be- 
come larger in ttifsssw- 



5 SUMMARY 

We have analyzed the complex, multimodal HB of NGC 2808 
as a prototype for a series of studies of hot HB stellar pop- 
ulations. We found that a combination of optical and far- 
UV magnitudes is the most powerful tool to efficiently an- 



Table 1. Mass (in solar units) at the MS-TO (M TO ) and mean 
value along the HB ((Mhb)), for each of the three HB sub- 
populations. AM represents an estimate of the mean mass lost 
during the RGB phase. (T e / Apeak ^ s * ne effective temperature of 
the peak of the distribution of each HB sub-population. 



POP 


Y 


(Mhb) 


M TO 


AM 


(-^cff )pcak 


RHB 
BHB 
EHB 


0.248 
0.30 
0.40 


0.84 
0.76 
0.627 


0.69 
0.565 
0.479 


0.15 
0.195 
0.148 


5,400 
12,600 
24,000 



alyze multimodal HBs with extremely blue components, in 
terms of both comparisons with theoretical models and di- 
rect star count analyses. By following this approach, the 
presence of He-enhanced populations along the blue part of 
the HB can be inferred without making assumptions about 
the RGB mass-loss. This cannot be done using only opti- 
cal CMDs because HB sequences with different initial Y 
overlap for T e g > 10, 000 K, where the HB becomes al- 
most vertical because of the large increase of the bolomet- 
ric corrections with T e ff. We have used here the (m,F336W, 
W.F336W — rnF555w) CMD to study the RHB and BHB stars 
and the {m f isa bw , rriFieoBW — ttifsssw) plane for hotter 
sub-populations (EHBs and BHks). BHB stars have been 
used as "reference" population to check the completeness 
level of the samples observed in the two different planes. 

From the comparison of theoretical ZAHB sequences 
with observations, we confirm that the peculiar HB of 
NGC 2808 cannot be reproduced by a population with a sin- 
gle initial He abundance, instead we need the combination of 
three different populations with different helium abundances 
ranging from Y = 0.248 to 0.40. Our models use Y = 0.248, 
0.30, 0.40 which are consistent with those inferred by Piotto 
et al. (2007) from the observed split of the cluster MS. Us- 
ing synthetic HB simulations we find that the CMD of the 
RHB sub-population is well matched by models with stan- 
dard initial Y values (Y = 0.248 in this case) and an average 
mass (Mhb) = O.69M0. The BHB stars can be reproduced 
by models with Y = 0.30 and (Mhb) = O.57M . The EHB 
sub-population is matched by models with with Y = 0.40, 
that might be the progeny of the blue-MS population. These 
latter stars arrive on the ZAHB with a very tiny H envelope, 
that results in surface temperatures T c ff > 20, 000 K, with 
a total mass (Mhb) = 0.48Mq. By comparing the average 
mass of the three sub-populations along the HB with those 
at the RGB tip for a 12 Gyr isochrones, [Fe/H] = -1.31 
and the appropriate values of Y , we find that the popula- 
tions have experienced an average mass-loss in the range 
(AM) ~ 0.17 ± 0.02 M Q . On the basis of HB evolutionary 
tracks and synthetic HB modeling we find that the HB of 
NGC 2808 is one of the most extended in temperature, with 
stars spanning the interval 5, 000 < T e g < 50, 000 K. These 
temperatures, estimated from stellar evolutionary models, 
are consistent with those derived from spectroscopic obser- 
vations by Moehler et al. (2004) , but significantly larger than 
those determined from optical CMDs. 

UV magnitudes for stars with T e g > 11, 500 K are af- 
fected by the levitation of heavy elements greatly increasing 
atmospheric opacity. We have crudely modeled the effects 
of levitation by employing bolometric corrections appropri- 
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ate for metallicities much higher than that of NGC 2808. 
The resulting models agree quite well with the observations 
especially for the BHB. Still, because of this approximate 
treatment our values for Y and the masses of the HB sub- 
populations are uncertain. However, our result that there 
are 3 populations with helium abundances spanning a range 
of roughly 0.15 in Y is quite robust. 

Even models with extremely large Y abundances are 
not able to reproduce the location of the BHk stars. This 
sub-population has been interpreted here as hot-flashers, 
the result of He-flash ignition along the He-core white dwarf 
cooling sequence. Why these stars are present in a few mas- 
sive clusters and absent in other clusters with comparably 
long blue tails is a challenge to stellar evolution theory. 

A comparison between HB and MS sub-populations 
relative ratios reveals that while (-/Vrhb /-^VEHB)obs is com- 
patible with (-/Vrhb /^Vehb)cxp, it is non-trivial to ex- 
plain the discrepancy between the values obtained for 
(A r RHB/A r BHB)obs and (iVRHB/-/VBHB)exp. The four sub- 
populations (RHB, BHB, EHB and BHk) show the same 
radial trend for r < 75", in agreement with previous findings 
(see BOO and 109). However, there is some evidence that the 
sub-population ratios might vary over larger radial distances 
(109). Studies of radial variation of the various populations 
covering the entire cluster coupled with detailed modeling 
of dynamics would be valuable. 

This paper is the first in a series which will make 
a detailed comparison between theoretical models and 
observations of blue tail HB clusters. 
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